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ABSTRACT 

This paper provided the first opportunity to quantify the angular stiffness and equivalent viscous 
damping coefficients of an axial spline coupling used in high-speed turbomachinery. A unique test 
methodology and data reduction procedures were developed. The bending moments and angular 
deflections transmitted across an axial spline coupling were measured while a non-rotating shaft 
was excited by an external shaker. A rotordynamics computer program was used to simulate the 
test conditions and to correlate the angular stiffness and damping coefficients. In addition, 
sensitivity analyses were performed to show that the accuracy of the dynamic coefficients do not 
rely on the accuracy of the data reduction procedures. 

INTRODUCTION 

Internal friction as a problem in rotor-bearing systems was first recognized in the mid- 1920s when 
manufacturers began shrinking disks onto supercritical rotors. Newkirk [1] observed that at speeds 
above die first critical speed, the rotor would enter into a violent whirling in which the rotor 
centerline processed at a rate equal to the first critical speed. Kimball [2] suggested that internal 
shaft friction could be responsible for the whirling. Since then, numerous contributors [3-6] have 
provided additional insight into the effects of internal friction. 

Today, with the increased power densities of advanced turbopumps and the resultant increased 
rotor flexibility, internal friction is seen to be a potential source of problems. Turbopumps such as 
die Space Shuttle Main Engine (SSME) High-Pressure Oxidizer Turbopump (HPOTP) are of 
built-up design with many joints, fits and areas for friction-induced excitation if slippage takes 
place. These rotors operate above flexible bending critical speeds and the forces on the rotors are 
very large, which tends to encourage joint slippage and friction force generation. The rotor 
instability resulting from this type of vibration is usually cured by adding special dampen, 
redesigning rotor fits, or modifying die bearings, which does not allow die development of a true 
understanding of the mechanisms of internal friction. In addition, the practicing design engineer 
finds very little readily available information on how to model the various couplings which make 
up the complete rotor system. Often, a rotordynamicist must make an intelligent guess on the 
dynamic properties of the coupling in order to perform the analysis. 

In the design charts of flanged and curvic couplings provided by Bannister [7], flexural stiffness 
values were provided but no damping data were available. Marmol et. al. [8] developed a 
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mathematical model to predict the lateral and angular stiffness and damping coefficients of spline 
couplings. Based on these coefficients, they calculated the natural frequencies of a rotor system 
and compared them to the experimental date. Recently, Walton et. al. [9] developed three 
nonlinear analytical models representing interference fits and axial and curvic spline couplings, and 
integrated them with a transient rotordynamics computer program analysis. Results from this work 
show the predictability of rotor system instability due to internal friction, which has been 
experimentally observed, especially for axial splines. 

It has been shown that measuring only the natural frequencies of a rotor system cannot accurately 
verify models predicting angular stiffness and damping coefficients of axial spline couplings [10]. 
However, no direct measurement has been reported, such as measuring the angular displacement 
and moment. This paper describe a unique test methodology and data reduction procedures that 
have been developed determining the angular stiffness and damping coefficients of an axial spline 
coupling. The bending moments and angular deflections transmitted across the axial spline 
coupling were measured while a non-rotating shaft was excited by an external shaker. The angular 
stiffness and equivalent viscous damping coefficient were determined by the interactive use of 
rotordynamics computer program simulations and experimental measurements. In addition, 
sensitivity analyses were performed to demonstrate that the accuracy of the data reduction 
procedures does not affect the accuracy of die results. 

EXPERIMENTAL TEST SETUP 


The test facility was built to simulate an SSME HPOTP prebumer spline. Figure 1 shows the 
overall test setup and Figure 2 is a functional block diagram of the test setup. The shaft consists 
of two offset male axial splines and an external female spline sleeve. The shaft was excited by an 
external dynamic force which was driven by an MTS hydraulic force control system in the X- 
direction (horizontal). An external damper was installed but was not active while the bending 
moment and angular deflection were measured. The damper was used to limit the shaft orbit while 
the system was traveling through its unstable natural frequencies. The shaft total length was 1 .524 
m, supported cm two deep-groove ball bearings 1.357 m apart. Total shaft weight was 155.4 N 
without the spline sleeve installed and 222.4 N with the sleeve installed. The shaft diameter 
outboard of the splines was 37.34 mm, while the center section between the two splines was 29.97 
mm in diameter and 57. 15 mm in length. The active length of each shaft external spline was 38. 1 
mm while the spline sleeve active length was 133.4 mm in order to span both of the shaft external 
splines. A close-up of the axial spline coupling is shown in Figure 3. Both shaft splines were 
identical but were machined 0° 18' (±15") offset from each other. The static torque preload that 
resulted was approximately 678 N-m. This torque value is representative of the SSME HPOTP 
prebumer spline. The two shaft splines and the mating internal spline sleeve were fabricated 
according to the following criteria: 


Type 

Number of Teeth 
Pressure Angle 
Pitch 

Pitch Diameter 
Base Diameter 
Surface Finish 
Coating 


Fillet root side fit 

42 

30° 

20/40 
53.3 mm 
46.1937 mm 
32v 
MoS, 
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Kulite ADP-350-300 semiconductor strain gages and specially fabricated high-resolution 
capacitance displacement probes were used to measure the bending moments and angular 
displacements, respectively. Four pairs of strain gages were installed mi the shaft surface at two 
axial locations (see Figure 3). At each location, one pair of strain gages was located at 0°/180° 
with respect to the key phase, and another pair was located at 90°/270°. Eight capacitance probes 
were installed to measure the angular deflection of the shaft and spline sleeve at four axial 
locations, where a disc was attached to each end of the spline sleeve and two discs were attached to 
the shaft outboard of the spline sleeve. At these four axial planes, axial motion was measured in 
line with the shaker stinger and at 90° to the stinger, and was converted to the angular deflections 
of that plane. In all cases, precise angular position of the instrumentation is required to ensure the 
integrity of the data. Circumferential alignment of the bending moment strain gage bridges was 
maintained to within 0. 125°. Likewise, foe angular deflection capacitance probes were positioned 
to the same accuracy from foe shaft centerline plane. In addition to foe bending moment and 
angular deflection instrumentation, shaft lateral motion in two orthogonal directions, input force, 
and acceleration of foe input force ram head were measured. While the testing discussed herein 
relates to nonrotating testing, the instrumentation setup is capable of being used as the shaft is 
rotated and excited nonsynchronously by the shaker [1 1]. 

As shown in Figure 2, foe data was acquired by a computer digital data acquisition system. For 
each test condition, foe input force, input acceleration, frequency, shaft motions, axial capacitance 
probe displacements, and bending moments were measured. To ensure a clear reference signal for 
all of foe data, the oscillator used to drive foe MTS hydraulic actuator was run through a 
synchronous tracking filter and subsequently converted to a square wave signal with a sharp 
leading edge. This square wave signal was then used as foe reference excitation timing signal for 
all data acquisition. 

Prior to conducting foe shaker testing, system qualification included a series of low-frequency tests 
where hysteresis loops of the input strain were plotted against foe bending moment beneath the 
spline sleeve [11]. Once the data instrumentation was qualified, a series of tests was run to acquire 
bending moment and spline angular deflection data. During the test, data was taken with foe shaft 
oriented so that the 0°/l 80° strain gage pair was directly in line with foe shaker stinger. Then, foe 
orientation of foe shaft was rotated so that the 90°/270° strain gage pair was in line with foe 
stinger, foe data was retaken under the same conditions. During foe data taking, foe off-axis 
(orthogonal to foe stinger plane) bending strain gage data was reviewed and the system alignment 
was adjusted by rotating the orientation of the shaft manually 1° or 2° to minimize the magnitude 
of foe signal in that plane. Once the out-of-plane bending strain signal was minimized, it was 
judged that the system was properly aligned. 

DATA REDUCTION PROCEDURES 

A multi-level rotordynamics model of foe test system was developed for a finite-element-based, 
steady-state linear rotordynamics computer program [12]. The main shaft was divided into 41 
elements, and each element was connected by two stations. The axial spline sleeve was modeled as 
the second level and had 13 elements. The connections between the main shaft and spline sleeve 
were modeled as two identical bearings. Figure 3 defines the station numbers in this region and 
shows the locations where foe instruments were installed. Each bearing was assumed to have non- 
zero lateral stiffness, angular stiffness and angular damping, but with zero lateral damping. To 
simplify the test data reduction procedures, all stiffness and damping were assumed isotropic, and 
cross-coupling terms were neglected. Introducing anisotropic terms would lead the dynamic 
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coefficients to be time dependent [10]. Lateral stiffness value was determined by a model 
developed earlier [9]. The rotordynamics model was tuned to match natural frequencies and mode 
shapes obtained during impact tests as closed as possible [10]. 

In the linear rotordynamics computer program, the change in the bending moments, AM, and 
angular deflections, A(j>, across the axial spline can be represented by the following equation in 
terms of the angular stiffness, K A , and equivalent viscous damping coefficient, B A [6]. 

AM,-K A M 1 +B A Ai 1> i-x,y (1) 

Physically, it is very difficult to measure the bending moments and angular deflections at the spline 
coupling locations (stations 25 and 61). Therefore, to approximate the bending moments and 
angular deflections between stations 25 and 61 (61-25), the original test procedure called for 
subtracting the bending moments at stations 23 from the bending moment at station 29 (29-23) to 
arrive at an approximation of the bending moment at the spline coupling. Similarly, the angular 
deflection at the spline coupling was determined by subtracting die angular deflection at station 30 
from that at station 65 (65-30). The initial approximated angular stiffness and damping 
coefficients were computed by Equation (1). Then, the bending moments and angular deflections 
at those stations were calculated by the rotordynamics computer program based on the 
approximated coefficients. An iterative procedure was developed to use both the rotordynamics 
computer program and Equation (1) interactively until results of the coefficients converged. 

Table 1 displays the calculated relative bending moments and angular deflections, both amplitudes 
and phases, across the spline coupling, stations 61-25, and at the locations to be measured, stations 
29-23 and 65-30, for a wide range of angular stiffness and damping coefficients as the shaft was 
shaking with 222 N at 50 Hz. We found that the coefficients were very sensitive to the accuracy of 
the phase of the relative angular deflection to be measured. The same phenomenon was found at a 
lower frequency. This is due to the fact that the shaft is much easier to bend near the outboard of 
the spline; therefore, the angular deflection of the shaft changes dramatically from station 25 to 
station 30 (see Figure 4). At the same time, the spline sleeve responded like a rigid body, and the 
angular deflections of the spline sleeve had almost the same values at both stations 61 and 65. 
Unfortunately, tire change of the shaft angular deflection compensated for the change of the angular 
deflection across the spline coupling caused by the damping, especially the phase. Therefore, an 
alternative data reduction procedure was required to determine tire dynamic coefficients. 

Table 2 shows the calculated absolute bending moments and angular deflections, both amplitudes 
and phases, at stations 23, 29, 30 and 65, for the same range of angular stiffness and damping 
coefficients as shown on Table 1, but with an applied load of 222 N at 5 Hz. We found that the 
absolute amplitude of the bending moment at station 23, the neck section between two splines, 
decreased as the angular stiffness increased, and the phase, with respect to the phase of a key 
location, of the bending moment at station 23 increased as the damping decreased. As the 
frequency increased, the amount of the change in both amplitude and phase due to the change in 
angular stiffness and damping increased. The amount of change was both measurable and 
distinguishable, providing an opportunity to use only the bending moment at station 23 to 
determine the angular stiffness and equivalent viscous damping coefficient. Other measured data 
act as checking points against which to compare the calculated data. 

The revised data reduction procedure is: (i) measure the bending moments at both stations 23 and 
29, and angular deflections at both stations 30 and 65; (ii) calculate these bending moments and 
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angular deflections by a rotordynamics computer program based on guessing values of K A and B A ; 
(iii) adjust K A and B A until the amplitude and phase of the bending moment at station 23 match the 
measured data; (iv) perform minor adjustment of K A and B A to minimize the sum of the relative 
errors due to the difference between the calculated and measured data at stations 29, 30 and 65. 

RESULTS AND DISCUSSION 

The tests were conducted at three levels of external shaking forces: 80.4, 162.6, and 221.1 N at 5 
Hz while the shaft was not rotated. At each force level, three data sets were taken, and we found 
that the variations of all measured data were within 1% in amplitude and 0.2° in phase. As the 
orientation of the shaft was rotated 90°, the same accuracy was achieved. 

Table 3 displays the calculated angular stiffness and damping coefficients of the three force levels 
using the revised data reduction procedures. The comparisons of the measured data and predicted 
values at those four stations are also shown in the same table. For all three force levels, the errors 
between measured data and predicted values are all within 5%. However, we expect that the errors 
may increase as the shaking force frequency increases. 

Note that for these three test cases, the angular stiffness decreases and damping increases as the 
force increases. This may be due to the fact that the relative motion between contact surfaces of 
the spline and sleeve increases as the shaking force increases. To determine the source and the 
mechanism of damping of the spline coupling, more studies are required, especially testing at 
different static preload torque levels of the spline coupling; however, this requires a permanent 
change of the hardware. 

The sensitivity analyses were performed for the case with 222 N shaking force at 5 Hz to study the 
effects of the accuracy of both the rotordynamics model as well as the measured geometrical 
configurations. On Table 4, which shows the results of the sensitivity analyses, it can be seen that 
the values of the various parameters were changed somewhat from the baseline values. The change 
occurred in ball bearing stiffness, spline lateral stiffness, K r ; the shaft diameters, D B and D 29 , 
at station 23 and 29, respectively; the length between two splines, L^,,; the installed bending 
strain gage locations at both station 23 and station 29; the orientation of the shaft, (j)^ and the 
misalignment of the external shaking force stinger. We found that all the parameters have little 
effect on the accuracy of predicting bending moments and angular deflections at stations 23, 29, 30 
and 65. The biggest error is caused by the misalignment of the orientation of the shaft; however, in 
our experiment, we were able to reduce the misalignment to within 2°. After the sensitivity 
analyses were performed, we were able to conclude that the accuracy of the angular stiffness and 
damping coefficient determination relied on the accuracy of the measured amplitudes and phases at 
all four stations, but was unaffected by the revised data reduction procedures. 

CONCLUSIONS 

A unique test methodology and data reduction procedures were developed in order to measure the 
bending moments and angular deflections transmitted across an axial spline coupling while the 
nonrotating shaft was excited by an external shaker. The angular stiffness and equivalent viscous 
damping coefficients were determined by the interactive use of rotordynamics computer program 
simulations and experimental measurements. Sensitivity analyses were performed to demonstrate 
that the proposed data reduction procedure has little effect on the accuracy of the results. 
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For the three shaking force levels tested at 5 Hz, the errors between the measured data and the 
predicted values are all within 5%. In addition, we found that for these three test cases the angular 
stiffness decreases and damping increases as the force increases. However, more tests are required 
to determine the source and the mechanism of damping of the spline coupling, such as the 
measurements for different static preload torque levels of foe spline coupling. 
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Figure 1 Experimental setup. 



Figure 2 Experimental test setup functional diagram. 
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Figure 3 Close-up of the axial spline coupling and instrumentation installed locations 
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Figure 4 Predicted angular displacements at 225 N external force 
(K a =1.13x 10 6 N-m/rad, B A =1.13xl0 4 N-m-s/rad) 
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Table 1 Predicted relative bending moments and angular deflections 
(222 N horizontal external shaking force at 50 Hz) 


JL 

B* 

EfSBHB 







Phase 


Phase 

warn 

Phase 

1 

Phase 

N-m/rad 

N-m-s/rad 

N-m 

mm 

N-m 

mm 

rad 

mm 

rad 

mm 

1.13E6 

1.13E2 

27.89 

■a 

■3S3 

0.5 

1.167E-4 


2.281E-5 

-1.3 


1.13E3 

28.36 

iNKXiM 

26.32 

4.4 

1.160E-4 

-1.6 

2.223E-5 

-13.0 


1.13E4 

34.94 

mm o 

33.78 

5.8 

1.046E-4 

-2.9 

0.907E-5 

-66.6 

1.13E7 

1.13E2 

35.03 

0.0 

33.85 

0.0 

1.040E-4 

■EO 


-0.2 


1.13E3 

35.03 

0.1 

33.85 

0.1 

1.040E-4 


2.995E-6 

-1.7 


1.13E4 

35.12 


mm* 

0.6 

1.039E-4 

-0.3 

2.868E-6 

-16.9 


1.13E5 

35.99 

0.5 

34.96 

0.6 

1.023E-4 

-0.3 

0.939E-6 

-71.7 

1.13E8 

1.13E2 

35.99 

0.0 

34.95 

HI 

1.023E-4 

0.0 

3.094E-7 

0.0 


1.13E3 

35.99 

0.0 

34.95 

0.0 

1.023E-4 

0.0 

3.094E-7 

-0.2 


1.13E4 

35.99 


34.95 

0.0 

1.023E-4 

0.0 

3.093E-7 

-1.8 


1.13E5 

35.99 

0.1 

34.96 

0.1 

1.023E-4 

0.0 

2.953E-7 

-18.4 


Table 2 Predicted absolute bending moments and angular deflections 
(222 N horizontal external shaking force at 5 Hz) 


K. 

b a 

My (23) 

My (29) 

4>y (30) 

♦y (65) 



Amp. 

Phase 

Amp. 

Phase 

Amp. 

Phase 

Amp. 

Phase 

N-m/rad 

N-m-s/rad 

N-m 


N-m 

msm 

rad 

1531 

rad 

E 31 

1.13E6 

1.13E2 

■S3 

-0.1 

51.45 

0.0 

2.434E-4 

HI 

6.097E-5 

Hi 


1.13E3 

12.78 

-0.2 

51.45 


2.434E-4 

-0.2 

6.098E-5 



1.13E4 

12.40 

-11.8 

51.44 

0.0 

2.41 IE-4 

-2.3 

6.102E-5 

Hnw 

1.13E7 

1.13E2 

3.51 

0.0 

51.34 

0.0 

2.077E-4 

0.0 

6.117E-5 

0.0 


1.13E3 

3.51 

-0.1 

51.34 

0.0 

2.077E-4 


6.117E-5 

0.0 


1.13E4 

3.51 

-0.6 

51.34 

0.0 


0.0 

6.117E-5 

0.0 


1.13E5 

3.42 

-5.9 

51.34 

0.0 

2.073E-4 

-0.4 

6.117E-5 

0.0 

1.13E8 

1.13E2 

2.38 

0.0 

51.32 

0.0 

2.033E-4 

0.0 

6.116E-5 

0.0 


1.13E3 

2.38 

0.0 

51.32 

0.0 

2.033E-4 


6.116E-5 

0.0 


1.13E4 

2.38 

0.0 

51.32 

0.0 

2.033E-4 

0.0 

6.116E-5 

0.0 


1.13E5 

2.38 

-0.1 

51.32 

0.0 

2.033E-4 

0.0 

6.116E-5 

0.0 


* With respect to 180° 
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Table 3 Test cases with horizontal external shaking force at 5 Hz 



Unit 

E9RH 

F=162.6 N* 

F=221.1 N* 

K a - calculated 

N-m/rad 

2.79 E6 

2.32 E6 

2.31 E6 

B a - calculated 

N-m-s/rad 

3.7 E3 

6.0 E3 

6.9 E3 

My(29) - measured 

N-m 

19.10 

38.61 

52.58 

My(29) - calculated 





M y ( 23) - measured 

N-m 

2.83 

6.06 

8.71 

My<23) - calculated 

N-m 




M x (29) - measured 

N-m 

0.39 (2.0%)@ 

0.75 (1.9%)@ 

1.59 (3.0%)@ 

Mj{(29) - calculated 


0 

0 

0 

M x (23) - measured 

N-m 

0.07 (2.5%)@ 

0.18 (3.0%)@ 

0.09 (1.1%)@ 

M x (23) - calculated 


0 

0 

0 

M y ( 23) phase 
• measured 


-1.43 

-2.6 

-3.2 

My( 23) phase 
•calculated 

degree' 

-1.4 

-2.6 

-3.2 

$ y ( 30) - measured 

pm 

4.72 

9.58 

13.18 

^(30) - calculated 

pm 


9.73 (1.6%) s 

|||||gggj 

4> y ( 65) - measured 

pm 

1.30 

2.59 

3.48 

<t» Y (65) • calculated 

pm 


2.67 (2.9%) s 



* Data arc average of 3 test cases, the deviations are: amplitude (± 1%), phase (±0.2°) 

# With respect to 1^(29) 

t 

Relative error with respect to measured data 
® Relative error with respect to measured data at M Y 
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Table 4 Sensitivity analyses of 222 N horizontal external shaking force at 5 Hz 



MJ23) 


M^) 

My(29) 

M x (29) 

<M30) 

<t> x (30) 

<M65) 

4 X (65) 


Amp. 

Phase 

Amp. 

Amp. 

Amp. 

Amp. 

Amp. 

Amp. 

Amp. 

K B ~1.75E7N/m 

0 . 01 % 

— 

— 

0.07% 

— 

0 . 8 % 

— 

3% 

— 

K R ±7.01E8N/m 

— 

— 

— 

— 

— 

— 

— 

0.3% 

— 

D 23 ±0.51 mm 

5% 

— 

— 

— 

— 

0 . 1 % 

— 

0 . 1 % 

— 


0.3% 

— 

— 

— 

— 

— 

— 

— 

— 

D 29 ±0.51 mm 

0 . 01 % 

— 

— 

0 . 02 % 

— 

1 % 

— 

3% 

— 

M 23 location ±2.5 mm 

1.3% 


— 

— 

— 

— 

— 

— 

— 

M 29 location ±2.5 mm 

— 

— 

— 

0.4% 

— 

— 

— 

— 

— 

fey ±5°* 

0.4% 

.... 

8.7% 

0.4% 

8.7% 

0.4% 

8.7% 

0.4% 

8.7% 

Force Misalignment 
±2.5 mm 

0 . 01 % 

HI 

— 

— 

— 

— 

§n§ 

— 

— 

Combination 

7% 

Bfil 

8 . 2 % 

0 . 1 % 

8 . 8 % 

0 . 6 % 

8 . 7 % 

6 % 

9.3% 


The base line case uses the following parameters: K A -2.03E6, B A »5.65E3, K B -1.23E8, K R «1.21E8, 
DJ3-29.97 mm, DJ9-37.34 mm, 1-22.24-57.15 mm, 4^-0° 

M J 

Experiment is able to reduce the misalignment of the orientation of the shaft to within 2° 
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